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Abstract 

Based on the Dirac approach we have developed the relativistic vision of the WKB method 
for centrally symmetrical potential with mixed Lorenz structure. We have obtained relativistic 
wavefunctions of light quark and the new rule of quantization containing the spin-orbit interaction. 
These gives us the possibility of finding the energy levels and decay width of hydrogen-like quark- 
systems. 
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1. Introduction 



The Dirac equation with potentials of a concrete type, for which possible to write the exact solution, 
O meet enough seldom. Most often for searching the solutions either numerical or asymptotic methods 
<—| >are used. In many theoretical and applied problems a possibility of obtaining the asymptotic solution 
P-i.allows to carry out the fullest analysis of a problem. Therefore hardly there is a necessity in detail to 
Q^explain importance of creating and investigating asymptotic methods of the solving the Dirac equation. 
The Wentzel-Kramers-Brillouin quasi-classical approximation (or WKB method) is one of basic and 
* most universal asymptotic methods of solving problems of theoretical and mathematical physics (see, for 
k> ■example, PQ EJ EJ E] ) , for which the exact solutions are either unknown or rather onerous. As is known 
5_j JHEIEIE] in case of a Coulomb field it has a high accuracy even for small values of quantum numbers. In 
5^ [contrast to the perturbation theory the given approach is not connected with a smallness of interaction 
and consequently has more wide applicability area allowing to study qualitative legitimacies in behaviour 
and properties of quantum mechanical systems. The new application of quasi-classical approach can 
be low-energy sector of QCD (the energy spectrum of hadrons, confinement study, decay widths of 
hadrons), where standard approaches based on perturbative theories are inapplicable in consequence of 
the fact that interaction between quarks in this area is not small. 



2. The WKB method for the Dirac equation in a centrally 
symmetrical field 

Obtain the formulae of the quasi-classical approximation for solutions of the Dirac equation with po- 
tential of centrally symmetry having mixed Lorenz structure. Corresponding to this we search wave- 
functions of the stationary states (in standard representation) in the form of 

x5,- r -i( F(r)tojim(p) \ m 
I iG(r)n jVm (n) i' [L) 



where Q - the spherical spinor, j and m - the total angular moment and projection of j (J = I ± 1/2), 
/ - the orbital moment (/ + /' = 2j), n = r/r. 

After separation the angular variable the Dirac equation with centrally symmetrical potential con- 
taining both vector V(r) and scalar S(r) parts, takes form (ft = c = 1) 

dF k 1 

— — I — F — — [(E — V(r)) + (m + S(r))} G = 0, 
dr r n 

^-- r G+\[{E- V{r)) - (m + 5(r))] F = 0, (2) 

where F and G are the radial functions, k = =F (j + 1/2) for states with j = / ± 1/2, i£ - the energy of 
level. 

For finding the quasi-classical solutions of the system (J2J) it is convenient to write equations (J2J) in 
the matrix form 

<f-±Z>* W "^ /r E-V(r) + m + S(r)\ 

ft ' \GJ' ^ -£ + V(r)+m + S'(r) ft£;/r j ' ^ 

Here a Planck constant ft is remained in the explicit view, the prime denotes the derivative with respect 
to r. The solution of the matrix equation © we shall look as the formal expansion in powers of ft: 

* = <pexp\Jydrj, y{r) = ~y-i(r) + y (r) + hy x {r) + h 2 y 2 (r) + <p(r) = f>V n) ( r ), ( 4 ) 

where the upper (lower) component ipffl ( ipo^) of the vector ip^ corresponds to the radial wavefunction 
F(G). Having substituted (HJ) into © and equated to zero the coefficients of each power of ft, we arrive 
at the recurrent system 

(D - y_0 = 0, 

n 

(D - y_0 y^ +1 ) = ^ (n) ' + E tt^V, n = 0, 1, . . . (5) 

fc=0 

From the first equation of system (0) follows that y is eigenvalues, (p n = (pi is one of (right) eigenvectors 
of the matrixes D. Eigenvalues y_j (r) = \ are roots of the secular equation det (D — y_ x ) = 0. 
Corresponding eigenvectors fi can be find in explicit form at diagonalization of the matrix D — Xf 



U2 

y_i = ±y (m + 5) 2 — (E — V) 2 + — = ±q = X u (6) 

/ m + S + £-y\ / A^-fcr- 1 \ 

^ = jfer^ + A, J=^m + S-£ + V J" (7) 

Here index i = ±, Ai(r) and A 2 (r) are some functions, which will be determined below. 

Since matrix D is not symmetrical, together with right eigenvectors <pi, it is necessary to enter left 
eigenvectors ft. They define by conditions 

ft (D - X { ) = 0, ft ^ <pf, (8) 



ft = Si ( m + S-E + V, kr" 1 + Xi ) = B 2 ( Xi - kr~\ m + S + E-V) 



(9) 



Furthermore, left and right eigenvectors the mutually orthogonal 



(<fh <Pj) = (&)«* (Vj) a = const • S i:j . (10) 



0=1 



Remind that (pi there is vector-column, and Sij is the Kronecker symbol. For determination y Q , we take 
ifo = (pi in the first equation of (0) and multiply both parts of this equation by (pi on the left. Then 
due to (jHJ), the left part is zero and we get the equation for y , from which 

»<r) = -£H (ID 
We select functions A 12 (r), _B 12 (r), in (J7J), (J9J) so that the condition 

(<Pu = <Pi) ( 12 ) 

is realized. Then 

r 

J y dr = In y?i)~ 1/2 ( 13 ) 

and as a result we obtain 

* = ^ ^)~ i/2 ^ y ^ rfr j > ( i4 ) 

that is similar usual expression for quasi-classical wavefunction in nonrelativistic quantum mechanics 

~p _1 / 2 exp f±z J pdr^j . (15) 

Using the first two equation of system (J5J) by means of left and right vectors technique we find 
the terms yi, yo and ip^\ Solving next equations of the system (j3J) by the similar procedure one 
can sequentially find the terms 2/2, 1/3,..., </^ 2 \ V >••■ m the expansions But formulae for them 
are rather cumbersome, therefore in applications ones usually restrict themselves to only first terms. 
Actually the reason of this is the fact |T], [2] that the expansions in powers of H (jlj) in the general 
case don't convergent and are asymptotic series, the finite number of terms of which gives the good 
approximation for the wavefunction, if a parameter of an expansion (the Planck constant h) is rather 
small. 

It is always easy to satisfy the condition (fT2"|) . By substituting the expression (jZJ), (JHJ), into (fT2*|) we 
come to the equation 

AiB[ - A! x Bx (m + S) V + (E - V) S' 



A x Bi ~~ q(q±kr~ 1 ) 



(16) 



whene 



F 
G 



2q(q±- 
\ r 



-1/2 



exp < ± J qdr+ 



1 j< m + S)V>HE : V)S> U m + S + E-V\ 

2 J q(q±kr- 1 ) \\ kr' 1 ± q J y 1 



(here and hereinafter h = 1). Using the second form of writing the eigenvectors ipi and (pi (with 
multiplier A 2 (r) and B 2 (r) in ((7J) and leads to the following result: 




2q[qT- 
\ r 



(ro + S)V'+(E- V) S' 
q (q =p fcr -1 ) 




We note that q coincides (to within multiplier % = \/— 1) with radial momentum of the relativistic quasi- 
classical particle. The sign +(— ) in (JTTj) and (fTHj) corresponds to the solution increasing (decreasing) 
with increasing r. For the decreasing solutions (the sign — ) it is necessary to use the formula (|17p when 
k < and formula (|18|) when k > 0; for increasing solution - on the contrary. The choice of the suitable 
form of writing the solution is defined by that to value Q = q+\k\/r should be positive in below-barrier 
region. 

The effective potential 

E S 2 - V 2 k 2 

' ro 2m 2mr 2 ^ 

corresponds to the Dirac system (J2J). So U (r, E) looks like a potential with a barrier. We consider the 
most general case, when the effective potential is barrier type. Then wavefunction has the different view 
in the three regions: 1) potential well r < r < r_ {q 2 < 0); 2) the below-barrier region r_ < r < r + 
(q 2 > 0); 3) the classically allowed region with continuum spectrum r > r + ( q 2 < 0); 4) classically 
forbidden region r < ro (q 2 > 0), where ro, r_ and r + are turning points. In following section we 
consider behaviour of the solutions in these regions. 

3. The wavefunction of the Dirac particle in classically allowed 
and forbidden regions 

The wavefunction of state has the different form in the various regions. 

I. The region ro < r < r_ is classically allowed; there the wavefunctions (|17|). (|18p oscillate 



E-V+m+S 
P 



nl/2 



cos©!, G = Cfsgnk 



E-V-m-S 
P 



nl/2 



COS 0- 



where 



p(r) = ^J(E-V) 2 -(m + S) 2 - (k/r) 2 
is the quasi-classical momentum for the radial motion of a particle, Cf - normalization constant. 



0x= / [p + ^dr + ^ 
V'-S' 1\ 



e, 



kw\ . 7T 

p H \dr + —. 

pr J 4 

V' + S' 



(20) 
(21) 

(22) 



w 



w 



+ 



1 



2\m + S + E-V rj 2\m + S-E + V rj 

Signs ± correspond to values k > and k < 0. If the level width F is small (it will be shown later) the 
wavefunction of quasi-stationary state can be normalized on a single particle localized in the region I, 
neglecting its penetrability into the classically forbidden regions r < r and r > r_ pQ 



F +G ) dr = 1 



(23) 



Here cos 2 Bj (r) can be replaced with average value 1/2: 



E-V{r) 
P i r ) 



-1/2 



dr 



2s 1/2 

Tj 



(24) 



where T- the frequency period of a relativistic particle inside a potential well. We note that in turning 
points r and r_ equation 



E — V 



(m + S) 2 + A;/r 2 



1/2 



is true and P — V > m + S* in region I. 

II. The below-barrier region r_ < r < r + is classically forbidden. Here p = iq, and quantities q, y_i 
and yo are real. As known [I] the wavefunction should exponentially decreases inside of this region. So 
the solutions of the Dirac system (J2J) in the below-barrier region 
for k < are 



cxp 



for fc > 



c 9 - 



exp 



-9 + 



jm + S) V' + (E-V) S r 

2^g , 



(m + 5 1 ) V" + (E - V) S n 
2^Q / 



dr 



-Q ) 

m+S-E+V h 



(25) 



dr 



m+S+E-V 

-Q 



(26) 



III. The result for region of continuum (r > r + ) is the most interesting. Here the wavefunction corre- 
sponds to the divergent wave (taking off particle): 
for k > 



Ci 



for k < 



exp 



cxp 



ip — 



r+ 



{m + S) V' + (E-V) S r 
2^P / 



(m + S) V + (E - V) S r 
2^P , 



dr 



iP 

m+S-E+V 



(27) 



dr 



m+S+E-V 
iP 



(2f 



where P = p + i|/c|/r. 

IV. In classically forbidden region (r < ro) the wavefunctions are of the form: 
for k < 



for k > 



- fc/r) 



^/g (q + k/r) 



exp 



2 + 



exp 



(m + g) V" + (E - V) S r 
2q (q - k/r) 



{m + S) V + (E-V) S r 
2q (q + k/r) 



dr 



dr 



E-V+m+S 
q — k/r 



q + k/r 
m+S-E+V 



(29) 



(30) 



The formulae (|20 p -()30j l include the whole range of values of r, except for neighborhoods of turning 
points. For bypass of these points and matching the solutions we use the usual method [2J. Closely to 
the r_ and r + the system (J2J) reduces to the Schrodinger equation with the effective potential linearly 
depending on r — r±, the solution of which expressed through the Airy function; one can match by the 
more elegant Zwaan method PQ, So the relation between the constants in various regions is of the 
form 



ct 



± 



E -V (r_) + m + S(r^ 
\k\rZ 1 



cxp 




2Ci 



\k\ rn 



E-V (r ) + m + S (r ) 



m + S) V' + (E-V)S' S 



2qQ 

-l) n Cf. 



dr 



sgnfc 

~ir~ 



(31) 



Let us find the equation determining the real part of the level energy E and width V of quasi- 
stationary levels E = E r — iT/2. Neglecting the barrier penetrability from ()23j) we obtain the quanti- 
zation condition: 



[ (p-\ | dr 

J \ pr I 

r x ' 



n + - ) 7T. 



(32) 



where n — 0, 1, 2, . . . is the radial quantum number. 
Pass to calculation of the level width 



-21m [G*F\ 



(33) 



Having used the explicit expressions ()27|) . (|2*Hj) for F and G functions, relation between normalization 
constants (}3*Tj) and expression for (j2U), we get the tunnelling probability: 



T 



cxp 




dr 



(34) 



Though the formulae ()23|) - (j3*4*|) essentially differ from the formulae of nonrelativistic quasi-classics (in 
particular, by relativistic expression for quasimomentum p, taking into account spin-orbit interaction, 
and additional pre-exponent multiplier) and more complicated from them. Their application to con- 
crete problems does not meet difficulties, since all quantities in functions F and G are expressed in 
quadratures. 

4. Approbation of results obtained 



4.1. For testing elaborated by us version of the WKB method, get the known in atomic physics the 
Bohr-Sommerfeld formula. We choose V(r) = — ^f- and S(r) = (a ~ 1/137 - the fine structure 
constant). Then the formulae (J21)) . (|22j) are of the form 



p(r) 



r 



.1/2 



m 



(k/r) 



1 

w = - 
2 



aZ/r 2 



E + ^+m 

r 



Having calculated the integral from ro to r_ 



EaZ =f J(EaZ) 2 - (m 2 - E 2 ) (k 2 - (aZ)' 



E 2 



in quantization condition 



we arrive at the expression 
EaZ 



a/ m 2 — E 2 



n + v k 



(azy 



whence we get expression, which is similar to the Bohr-Sommerfeld formula 

/ aZ \ 

E = ±m 



1 + 



n+^Jk 2 - (aZf 



4.2. One more rather interesting case, in which it is possible to get exact solution of the Dirac 
equation with potential of the type V(r) = — ^ S(r) = — — , exists. In the same way that the Coulomb 
potential V(r) is derived from the exchange of massless photon between the nucleus and the leptons 
orbiting around it, the scalar potential S(r) is created by the exchange of massless scalar mesons. The a 
meson frequently quoted in the literature has a very high mass and therefore the corresponding potential 
has a very short range. In this case the formulae (|2ip. (|22|) take the form 



a 



E+-) - [m (k/r 



p(r) = 

Quantization condition (|32|) gives following result 

Eol + ma! 



1/2 



1 

W= 2 



(a — a')/r 2 



E + 



+ m r 



\Jm 2 — E 2 

there was the integration from r$ to r_ 



n + 7, 7 = ±Vk 2 — a 2 + a 



r> 



r o,- 



Ea + ma' =F J (Ea + ma') 2 — (m 2 — E 2 ) 7 2 



m^ 



E 2 



Thus we find the formula coinciding with the expression obtained in 



E = m 



—aa 



a 2 + (n + 7)' 



aa 



s a 2 + (n + 7)' 



a' 2 — (n + 7)' 
a 2 + (n + 7) 2 



1/2) 



> . 



4.3. We consider the potential V (r) = S (r) = ar 2 /4, a is the force of an oscillator. In this case the 
formulae (|2~T|) . (|2*2*j) is given by 



p(r) 



E 2 -m 2 -(E + m) % 2 - (k/rf 



1/2 



w 



1 

2r' 



We integrate in quantization condition (32) from r to r_ 



2aA; 2 
+ m 



1/2 



and gives following result 



(E-m) 



lE + m ... 1 

\k H — sgnA; = 2n, 

2a 1 1 2 6 



for k < get 

(2 \k\ + l)y/a-(E-m) \j2 (E + m) = -Any/a. 
The last equation for energy is cubic, its real solution is 

E nM = (2m + 8 ■ 2 2 / 3 m 2 A~ 1 / 3 + 2 2 ^A 1 ^) /6, 

where A = -B + y/B 2 - 1024m 6 , B = 32m 3 - 27a (1 + 2k + An). This expression completely coincides 
with result obtained in jjj. 

5. Description of the energy spectrum in two-quarks systems 
with Cornell potential 

Relativistic description of the bounded states always was an important problem of nuclear physics and 
elementary particle physics, in which relativistic characteristics of the light quarks play the important 
role. Our world mainly consists of the light u and d quarks (the protons - uud and neutrons - udd). 

Quantum chromodynamics (QCD) based on principles of the quantum field theory is justly consid- 
ered as the most consequent approach to solving this problem. However, standard perturbative QCD 
gives rather reliable recipes of the calculation of various characteristics only for description so named 
"hard" processes which are characterized by large transmitted momentum, and not applicable for cal- 
culation of the characteristics which are defined by "soft" processes (the mass spectrum, confinement of 
quarks, decay widths of hadrons). In the same time nonperturbative effects basically define the nature 
of forming bound states of interacting particles. Confinement is the result of circumstance that, unlike 
quantum electrodynamics in which interaction mediators - photons - are electroneutral, exchange par- 
ticles in quantum chromodynamics - gluons - possess non-zero colour charge and therefore can interact 
one with other. Thus confinement is not nested in the framework of the perturbation theory. On ac- 
count of this in the given time from principles of quantum chromodynamics the structure of interquark 
forces cannot completely be defined. 

Just study of relativistic properties of quarks systems, namely of effects caused by spins of quarks, 
enables to improve both the form of confinement part of potential, and Lorenz structure of potential. 

In the work 8\ starting from the QCD Lagrangian and taking into account both perturbative 
and nonperturbative effects, using the method of vacuum correlators the Dirac equation is obtained 
(rigorously for the Coulomb interaction and heuristically for the confining potential) for a system 
consisting of a heavy antiquark (quark) Q (b and c quarks) with the mass ni2 and light quark (antiquark) 
q (u, d, s - quarks) with mass much less m\. 

Such quark systems (B + = bu, B° = bd, B° s = bs, D° = cu, D~ = cd, D~ = cs) are QCD analogues 
of the hydrogen atom and thus its study is of fundamental importance. From the theoretical viewpoint 
the interest in heavy-light systems stems from several considerations. 

First, in the limit of one infinitely heavy quark, one hopes to get the dynamics of a light quark in 
the external field of a heavy one. That would be similar to the picture of the hydrogen atom. 

Second, since the external field is time-independent, one may hope to obtain a static potential in 
QCD together with spin-dependent forces, as has been done for heavy quarkonia. 

Third, in the heavy-light system one may study how the chiral symmetry breaking (CSB) affects 
the spectrum. When one quark mass is vanishing, in the chiral symmetric case the spectrum would 
consist of parity doublets, and the CSB would lift the degeneracy. 

Fourth, using the Dirac equation we implement explicitly relativistic dynamics and can study rela- 
tivistic properties of the spectrum, e.g., in the case of a vanishing quark mass, and also relativistic spin 
splitting in the spectrum. 



In the work [8] it is explicitly shown that a reasonable spectrum in Dirac equation occurs only when 
the confining potential is a Lorentz scalar and the Coulomb potential is the fourth component of a 
4-vector. In this case the scalar potential breaks explicitly chiral symmetry, and states with opposite 
parities are not degenerate. In the case of a vector confining part, only quasi-stationary states are 
found. 

We consider the most general type of the confining part, which taking into account also cases 
considered in jH]. Let us suggest that the static quark potential has vector and scalar properties of 
Lorenz transformation: 

V NR (r) = V v (r)+V a (r). (35) 

The quantity V means that the potential is a 4-th component of operator p^, S means that the potential 
has scalar nature. 

Following many authors EH] we assume the mixture of vector-scalar quark interaction potential 

V (r) = Voge + eV conU V s (r) = (1 - s)V conh (36) 

where Voge is the one-gluon exchange (OGE) term, Konf is the confinement part of the potential, e is 
mixing constant. Here the Lorenz nature of the one-gluon part of potential and the confining potential 
is different, the one-gluon potential is totally vector, while the confinement is a vector-scalar mixture. 
Very interesting reviews were done in fH\ IT2] concerning the choice of interaction potential. 

To simplify the calculations we consider the simplest case when e = 0. This corresponds to the 
one-gluon potential being totally vector and the confinement potential being totally scalar. 

So we choose the scalar and vector parts of the potential (|3l))l in the form EH HH U2] 

Voce (r) = --, V cori f (r) = /ir (37) 
r 

where (3 = /i = 0.18 GeV 2 , a s is the constant of the strong coupling with asymptotic freedom: 

12-k 1 

a s r = x, 38 

sK J 33-2N ln(Q/A) 2 V ; 

where Q is the transmitted momentum, A is taken to be equal to A = 0.14 GeV, N is the quarks 
flavours number. 

In work the WKB method was used for potential ()37|) in the Schrodinger equation. For solving 
of the Dirac equation system (j2J) with potential (f3Tj) . we shall use early elaborated by us version of the 
WKB approximation. 

Due to the confinement of quarks we are interested in only classically allowed region (r < r < r_, 
q 2 < 0, a potential well) that corresponds to only the discrete energy spectrum of the quark-antiquark 
system. Thus the quasi-classical momentum is of the form {h = c = 1): 



P ( r ) 







2 



E + — 1 — (m + fir) 2 — (k/r) 2 
r 



1/2 



(39) 



We represent left part of the quantization condition (|32j) in the form of sum of two integrals I\ and I2 

a a 

f P nil] 

h=J p(r)dr, I 2 = J —dr. (40) 



The integration is between the two classical turning points ro = b and r_ = a, which are real and 
positive roots (a > b > 0) of the equation 



r H r 



2m o £ 2 -m 2 2E(3 k 2 - (3 



-r — 



r + 



0. 



(41) 



Two roots of this equation are real and negatives (d < c < 0). Formulae ()4()j) can now be re-expressed 
in terms of a, 6, c and d as 



Ji = -/i 



a ( 3 , 2m ? 
/ V M 



2£/3 _|_ fc 2 -ff 2 
/i 2 /^ 2 r 



[(a — r) (r — 6) (r — c) (r — <i)] 1//2 



(42) 



and 



H(A-B) 



A + 



E + m 



dr 



2/i ) J (r-A) [(a - r) (r - 6) (r - c) (r - d)] 1/2 



- 5 



E + m 



dr 



2fi J J (r-B)[(a-r)(r-b)(r-c)(r-d)} 1/2 



(43) 



where A and which are roots of the quadratic equation [ir 2 + (E + m) r + /3 = 0. These equations can 
be represented in terms of complete elliptic integrals. Then quatization condition ()32|) gives following 
result 



y / '(a — c) (6 — d) 



(6 - cr 



+ jv 4 -n(-i/,x 



+ - 



(1 - v) (x 2 - u) 
b-c 



[N, ■ F (x) + N 2 ■ E ( X ) + iV 3 ■ n (z/, X ) + 



(jv 5 ■ n (i/i, x) - n 6 ■ n (z/ 1; x )) + n 7 ■ f ( x ) 



n + - i 7T. 



(44) 



where F (x), (x) and II (z/, x) are complete elliptic integrals of the first, second and third kinds and 
^ = jxV(6-c) _ 3^-2,-2^ + 3x^-c) _ ^ + 3^ ^ ^ + 



(b-c) 



C/jL 



No = -v 



3 ^3^(z/ 2 -2z/-2z/x 2 + 3x 2 )(6-c) 
m 2 CM 8 



(1 - i/) (x 2 - i/) 



3 (ft _ C ) (z/ 2 _ 2zy _ + 3x 2 } 2 ^ (1 _ „) (x 2 _ y) / + ^ _ ^_ m 2 - 



(1 - u) (x 2 " I/) 



(6-c) 



+ [m + -cfi) (v 2 - 2i/ - 2z/ X 2 + 3x 2 ) - f (6 - c) (3 X 2 - (x 2 + l) v) 



(1-u) ( X 2 -u) (k 2 -(3 2 ) 



(b-c) 



bcfi 



X 



(a — b) (c — d) 
\(a-c)(b-d) 



a — b 



a — c 



A-c 
A-b 



B-c 
B—b^ 



Equation (|44j) is an implicit relation for the energy E n ^. 

Accordingly wave functions of light quarks in classically allowed region are defined by formulae (J2*U|) . 
where normalization constant is of the form 



li J (a — c) (b — d) 



2[{Ec + (3)-F{ X )+E{b-c)-n{u,x)\ 



1/2 



(45) 



and phases 61 and 02 are 



6i 



y / \a — c) {b — d) 



(b-cf 



(1 - u) (x 2 ~ u) 

c 



N 1 -F(<j ) ,x) + N 2 -E( ( f ) ,x)+N 3 -Tl( ( j ) , u, X ) + N 4 ■ U ( 0, -u, X ) + L 



+ 



k 

+- 



b — c 
A-B 



(N 5 ■U(<f>,u 1 ,x)-N 6 -U (0, v 2 , x)) + N 7 ■ F (0, X ) 



71 

+ 4' 



(46) 



e 2 



y / '(a — c) (6 — 0?) 



(6-c) 



(1 - 1/) (x 2 - 1/) 



■ x 



iv 1 -F(0, x ) + iv 2 . J E(0,x) + iV3.n(0,z/,x) + iv 4 .n^,^,x) +l 



fc 



6-c 
A 1 -B i 



(N s ■ n (0, v{, x) - N 9 ■ n (0, z/*, X )) + iV 10 • F (0, X ) 
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Here 
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6. Quasi-independent quark model of hadrons 



One of simple of meson quarks models, which, despite of singleness, gives the good agreement with 
experimental data, is the quasi-independent quark model [T4*l ITS] . Instead of to consider a problem on 
the bound states of system of two quarks, intensive interacting one another, we shall consider quarks 
in a meson as independent particles, which move in mean or self-consistent field. 

In such scheme not only wavefunction of all system of two quarks has a sense, but also individual one- 
particle wavefunctions of each quark. These one-particle wavefunctions are described by the usual Dirac 
equation for particle in some field. Thus, quarks, which are included in a structure of a meson, occupy 
one-particle energy levels. To simplify calculations it is assumed that this mean field is spherically 
symmetric and quark motion in space is determined by motion of its center. Each of constituents 
interacting with the spherical mean field gets the state with a definite value of its energy. On the 
phenomeno logical ground the mass formula for J PC meson, which is n 2S+1 Lj state of gg-system, can 
be represented in the following form: 



M = ypl + m l + ypl + ™%. (48) 

In a system of the centre of mass (|pi| = [_p2 1 = p) after simple mathematical transformation we have 
obtained the relation 

^± 2 =ap + Pm 1 = E, (49) 

from which 



M = E + yE 2 -m\ + m|. (50) 
Here E is the eigenvalue of the one-particle Dirac equation (J2J). 

7. Results and conclusions 



The Dirac equation gives the spin-orbit splitting only into two levels while the experiment gives three 
levels. It is due to the fact that the individual quark spins can be composed in or 1 and the total 
angular moment can be 0, 1 or 2 in the case 1 = 1. But in the limit of an infinitely heavy quark mass 
m 2 — > oo degrees of freedom of quarks system are determined by the quantum states of the light quark 
with the total angular moment j = L + s. This gives two sets of levels j = 1/2, j = 3/2. For example, 
in the hydrogen atom, we don't take into account the nuclear spin: it is entered only as a hyperfine 
effect. For comparing our data with other data of therms with total moment we hold the following rule: 
P1/2 was considered as being an averaged mixture of l P\ and 3 Po, while P3/2 was considered as averaged 
mixture of 3 Pi and 3 P2 masses. This problem can have different solutions and it will be discussed in 
our next papers. 

Everyone dealing with the Dirac equation encounters the problem of Lorenz structure of the po- 
tential. Our calculations demonstrate that considered one-gluon potential being totally vector and the 
confinement potential being totally scalar lead to incorrect spin-orbit splitting of P-levels: M(P 3 / 2 ) < 
M(Pi/ 2 ). For correct description of sign of spin-orbit splitting in the Dirac equation it is necessary to 
take a part of confinement potential as vector that corresponds to (jHSj) - Hereinafter we will consider 
the case e 7^ 0. This conclusion concerning the Lorenz character of the confinement agrees with other 
authors [Td1[T7]. 
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